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Abstract 

We use the more intuitive approach due to Kolmogorov (and subsequently, Landau in his text on 
fluid dynamics) to calculate some third-order structure functions for quasi-geostrophic turbulence 
for the forward cascade of pseudo-potential enstrophy and the inverse energy cascade in quasi- 
geostrophic turbulence. 
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I. INTRODUCTION 



Quasi-geostrophic (QG) turbulence is a rather more realistic class of turbulent flow than 
the isotropic homogeneous 3D turbulence. It can be seen in the large scale flows on oceans 
and atmosphere; thus having profound geophysical and astrophysical significance. In the 
incompressible isotropic homogeneous 3D turbulence, Kolmogorov's four-fifths law[l] is a 
landmark in the theory of turbulence because it is a rare exact result. In three spatial 
dimensions, this law says that the third order velocity correlation function behaves as: 



where £ is the rate per unit mass at which energy is being transferred through the inertial 
range. The inertial range is the intermediate spatial region postulated by Kolmogorov where 
the large scale disturbances (flow maintaining mechanisms) and the molecular scale viscous 
dissipation play no part. This result is of such central significance that attempts are regularly 
made to understand it afresh and to extend it in other situations involving turbulence. If we 
consider the 2D turbulence [2], then in the inviscid limit, we have two conserved quantities 

- energy and enstrophy. This gives rise to two fluxes with the enstrophy flux occurring 
from the larger to the smaller spatial scales. The energy flux goes in the reverse direction. 
QG turbulence [3] stands somewhere in between 2D and 3D turbulences. In the inviscid 
limit, besides total energy, QG flows enjoy the possession of yet another conserved quantity 
which is conserved at the horizontal projection of the particle motion. We shall call this 
pseudo-potential vorticity to distinguish it from the potential vorticity that is conserved at 
a particle in a homentropic fluid. Deflning pseudo-potential enstrophy as half the square of 
the pseudo-potential vorticity, one would say that like 2D turbulence there arc two cascades 

— forward cascade of pseudo-potential vorticity and inverse cascade of energy — in QG 
turbulence which, however, is inherently three dimensional in nature. 

Recently, a paper[4] has calculated some structure functions in QG turbulence and has 
made illuminating revelation that isotropy in the sense of Charney[3] is useless in deriving 
the structure functions for QG turbulence. It has gone on to show that formulation of 

QG turbulence under the constraint of axisymmetry is productive. However, it criticized 
(though somewhat rightly) the ineffectiveness of use of tensorial quantities in the case of 
QG turbulence in deriving the results. Now, manipulating the tensorial quantities are at 
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the heart of the derivation of many important two-point velocity correlation functions and 
other ones [5]. The technique is very intuitive and straightforward. It has, recently, also been 
thoroughly used to find out various correlation functions for 2D turbulence [6] . In this paper, 
we shall closely (and trickily) follow the original Kolmogorov method put forward in details 
in the fluid dynamics text due to Landau and Lifshitz[5]; and repeated in the ref.-([6]), to 
derive structure functions in QG turbulence. The method has the extra advantage to being 
able to probe into the form for the two-point third order velocity correlation function in 
the forward pseudo-potential enstrophy cascade regime — this has remained untouched in 
ref.-([4]). 



II. TWO POINT THIRD ORDER MIXED CORRELATION FUNCTION 

First of all we shall briefly introduce the necessary equations (see ref.-([7]) for details). 
Let u be the three dimensional velocity of the fluid in a frame rotating with constant angular 
velocity Vt. The fluid body (such as ocean) is assumed to be of uniform density with free 
surface &X, z — ^{x,y,t). Suppose the bottom z — —H{x,y) is rigid. The shallow-water 
equations, then, are: 
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Here, h{x,y,t) = C{x,y,t) + H{x,y), v = [u^^Uy], V = {d,^,dy), = v.V and 
V = v{x,y,t). f is Coriohs parameter that is Taylor-expanded to write / = /o + f^y- Using 
the equations (2) and (3), one gets the relation: 
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Let us assume: a) Rossby number Ro <^ 1, b) Fractional changes in h are small and c) 
PL/fo ^ 1 where L is the horizontal scale of the flow. Imposing these three assumptions 
on the shallow-water equations one can modify the relation (4) to yield 



where q = V'^ip + f — f^ip/gHQ + Jq^Hq — H)/Hq {ip being g^/ fo) may be called pseudo- 
potential vorticity. Under the same assumptions, for QG flow, one also has the condition: 

V.v = Q (6) 

Now, the trick is to select an arbitrary two-dimensional plane in the QG turbulent flow such 
that the plane's normal is parallel to the vertical {i.e., along /) and impose the property 
of homogeneity and isotropy in the plane only. By the way, one must keep in mind that 
the so-called fundamental scale of 3D turbulence has its analogy as the horizontal length 
scale L for the case of QG turbulence and; the correlation functions to be derived for the 
forward cascade in this paper are valid in the range (which we shall call inertial range) that 
is much smaller than L but quite larger than the scale at which the dissipation is effective 
and, the structure function to be derived for the inverse cascade is valid in the range whose 
scale is larger than the scale at which energy is being fed in. As we shall consider fluid 
bodies of uniform density only, we shall take density to be unity and let p take over the 
task of representing position vector in the 2D plane. Yet another convention: The Greek 
subscripts used herein can take two values p and ± which respectively mean the component 
along the radial vector p and the component in the transverse direction. When we shall use 
the Latin subscript {e.g., a), it should mean that it can take one more value apart from the 
ones mentioned above: the third value '2;' would signify the vertical direction. Einstein's 
summation convention will be used extensively. Also, 

P = P2-P1, p"a = Pa/\p\, P;=l, p1 = (7) 

Now, if Vi and V2 represent the horizontal fluid velocities at the two neighbouring points at 
pi and p2 respectively then with similar meaning for qi and q2, one may define: 

K ^ (qm) (8) 
and, Q = {{q2 - qi){q2 - Qi)) (9) 

The angular brackets denote an averaging procedure which averages over all possible posi- 
tions of points 1 and 2 at a given instant of time and a given separation. Due to homogeneity, 
Q may be re-expressed as: 

Q = 2(g2) - 2K (10) 



For simplicity, we shall take a rather ideahsed situation of QG turbulence which is homo- 
geneous and isotropic on every scale in the plane. For the unforced case, the component 
of the correlation tensor will obviously be dependent on time, a fact which won't be shown 
exphcitly in what follows. As the features of local QG turbulence should be independent of 
averaged flow, the result derived below is applicable also to the local turbulence in the plane 
at scale p much smaller than the fundamental scale. 

Again, we define a two-point third order mixed correlation tensor in inertial range: 

Qs = {{V2S - Vis){q2 - qi){q2 - qi)) (11) 

Qs = 4Ks + 2Ms (12) 

where, 

Ks = {visqiq2) (13) 

and, Ms = {qiqiV2s) (14) 

Let due to isotropy and homogeneity, we can write following form for M^: 

Ms = M(p)p^ (15) 
d 

—Ms = {qiqid2sV2s) = (16) 

P2S 

^ Sm(p, + ^ = (17) 
dp p 

, constant „ 
Mip) = = (18 

P 

In the relation (16), we are using the expression (6) and in writing the relation (18) we have 
taken into account the fact that Ms should remain finite when p = 0. Relations (15) and 
(18) imply that: 

M5 = (19) 

using which in the relation (12), we get: 

Qs = ^Ks (20) 

From the equation (5), we may write for the points 1 and 2 respectively: 

d 

—q^ = -vi^di^qi (21) 
d 

—q2 = -V2^d2'yq2 (22) 



Multiplying equations (21) and (22) by q2 and qi respectively and averaging subsequently 
after adding, we get the following differential equation: 

= 2dsKs (23) 

where we have used the fact ds — —du — 825. Using relations (10) and (20) in the equation 
(23), one gets for the inertial range for the pseudo-potential enstrophy cascade in homoge- 
neous and isotropic QG turbulence (forced at an intermediate scale or unforced) in inviscid 
limit the following differential equation: 

^ Qp = -2e^p (25) 
In getting relation (25) from the equation (24), we have assumed the facts: 

1. \-§i{q^) = —£q, i-G., there exists a pseudo-potential enstrophy sink at small scales due 
to some dissipative force such as viscosity and Eg is the finite and constant dissipation 
rate of the mean pseudo-potential enstrophy. 

2. ~ due to quasi-stationarity. It may be supposed that the value of Q varies 
considerably with time only over an interval corresponding to the fundamental scale 
of turbulence and in relation to local turbulence the unperturbed flow may be re- 
garded as steady which mean that for local turbulence one can afford to neglect -^Q 
in comparison with the pseudo-potential enstrophy dissipation rate Eg. 

III. TWO POINT THIRD ORDER VELOCITY CORRELATION FUNCTION 

Again, one may define a rank two correlation tensor: 

BaS = {{V2a " Via){v2S " Vu)) (26) 

Isotropy and homogeneity in the plane suggests following general form for B^s 

B^s = A,{p)5^s + A2{p)p°ys (27) 

where Ai and A2 are functions of time and p. Making use of the relations (7) in the equation 
(27), one gets: 

BaS = S±± {Sas - pIpD + B,M (28) 



One may break the relation (26) as 



BaS = (viaVls) + {v2aV2s) " {viaV2s) " {v2aVl&) (29) 

Defining 

Ks = {viaV2s) (30) 

one may proceed, keeping in mind the isotropy and the homogeneity, to write 

BaS = {v^)SaS - 2ba5 (31) 

Again, having used the condition (6), one may write: 

dsBas = 

B^^^pB'^^ + Bpp (32) 

where the equation (28) has been used and prime (') denotes derivative w.r.t. p. Let's give 
yet another definition: 

Invoking homogeneity and isotropy in the plane once again along with the symmetry in the 
first pair of indices, one may write the most general form of the third rank Cartesian tensor 
for this case as 

bas,^ = C{p)S^sp; + D{p){5,sPl + ^a,pD + F{p)PaPsP"^ (33) 

where, C, D and F are functions of p. Yet again, the expression (6) dictates: 

_d_. _d_ _ 

C, 2D ^ 2D' 2D F' F ^ 

6a6^ OaS H Oo,S H ^-PaPi ^PaPS + -^PaPS + -^PaPS = (34) 

p p p^ p"* p^ p'' 

Putting q; = 5 in equation (34) one gets: 



where, it as been imposed that 60,5,7 should remain finite for p = 0. Again, using equation 
(34), putting a ^ S and manipulating a bit one gets: 

D^-^{pC' + C) (36) 

using which in relation (35), one arrives at the following expression for F: 

F^pC -C (37) 

Defining 

— 2(6a<5,7 + ^lS,a + ^07,5) (38) 

and putting relations (36) and (37) in the equation (38) and using relation (33), one gets: 

B^s^ = -2pC"((5o5p; + 5^spI + 5a^p1) + 6(pC" - C)p>^p° (39) 
^ Bppp^-6C (40) 

which along with relations (36), (37) and (33) yields the following expression: 



bas,y = -^S^sp; + ^{PB',,, + B,,,){5,5Pl + ^^..pI) " \{pBl, - B,,,)py°sP; (41) 



The equation (3) suggests: 

d 

—Via = -Vl^di^Via + fla€aa-,Vlj - gdia^l (42) 

d 

—V2S = -V2-yd2^V28 + ha^a6-yV2-y " gd28^2 (43) 

multiplying equations (42) and (43) with V25 and Via respectively and adding subsequently, 
one gets the following: 

d 

-Q^{viaV2s) = -dij{vi-yViaV2s) " d2j{v2jViaV2s) 
+^aaj{flaVljV2s) + ^a&-t{f2aV2'yVia) 

-gdla{ilV2s) - gd2s{^2Vla) (44) 

Due to isotropy, the correlation function (^1^/2) should have the form f{p)p/\p\. But since, 
da{CiV2a) = due to solenoidal velocity field f{p)p/\p\ must have the form constant x (p/|p|^) 



that in turn must vanish to keep correlation functions finite even at p = 0. Thus, equation 
(44) can be written as: 

d 

'O^^ctS = dj{ba'y,S + bs-y^a) + fo^za-yb-fS + fo^zS-fba-y (45) 

— * 

Here we have used the approximation: / = foz. Using equations (31) and (41), one can 
rewrite equation (45) as: 

1(9, 9, 1<9„ 1 d , . 

Note that the terms containing the Levi-Civita symbol vanish of the joint effect of the 
expressions (28) and (31); and the antisymmetry property of Levi-Civita symbol. As we 
are interested in the pseudo-potential enstrophy cascade, the first term in the L.H.S. is 
zero because of energy remains conserved in QG turbulence in the inviscid limit; it cannot 
be dissipated at smaller scales. Also, as we are interested in the forward cascade which is 
dominated by pseudo-potential enstrophy cascade, on the dimensional grounds in the inertial 
range Bpp (if it is assumed to depend only on eg and p) may be written as: 

^S,, = Fsgp' (47) 

where F is a numerical proportionality constant. Hence, using the relation (47), the equation 
(46) reduces to the following differential equation: 

which when solved using relation (40) imposing finiteness of Bppp for p — 0, one gets 

B,,, = -^P' (49) 

The relation (49) is the expression for the two-point third order correlation function in the 
isotropic and homogeneous plane of QG turbulence (forced or unforced) in the range of the 
forward cascade where there is no overlapping with energy cascade. Since F has not been 
determined one must confess that the equation (49) is just a scaling law at this stage. 
Now suppose the fiuid body is being forced at small scales i.e., energy is being supplied 
and the mean rate of injection of energy per unit mass is denoted by Su (assumed finite 
and constant). Let us focus on the inverse energy cascade. Then technically we have to 
proceed just as before to finally arrive at the differential equation (46). One obviously would 



set ^■§i{v'^) = invoking the hypothesis[3] that there should be equipartition of energy 
between potential energy and the energy content in each of the two horizontal velocity 
components in the plane. Lets also assume that §iBpp ~ in the inverse cascade regime 
supposing the forced QG turbulence to be in the state of quasi-stationarity. So we are left 
with the following differential equation: 

6^1 (p'BpJ = (50) 

=^ Bppp = +euP (51) 

where in the last step the integration constant has been set to zero to prevent Bppp from 
blowing up at p = 0. The expression (51) is the expression for the two-point third order 
correlation function in the isotropic and homogeneous plane of forced QG turbulence for the 
inverse energy cascade. 



IV. CONCLUSIONS AND DISCUSSIONS 



It has, thus, been again showcased how handy and useful the Kolmogorov- Landau ap- 
proach can prove to be. The results (25) and (51) naturally agree with what has been arrived 
at by Lindborg[4] earlier. Within the domain of the approximations made these results are 
exact, something worth getting as the literature of turbulence is comparatively barren as far 
as exact relations are concerned. However, the hypothesis of the equipartition of energy used 
in equation (50) is as questionable as the assumption of isotropy in the sense of Charney. 
This hypothesis needs to be put on more firm basis. In the closing, we hope that validity of 
the results derived will be checked both numerically and experimentally in near future to see 
if the approximations made for the homogeneous QG turbulence in this paper arc correct 
or not. Also, the fact that the structure functions for the inherently three-dimensional 
QG turbulence are more like that of the 2D turbulence than that of the 3D turbulence 
speaks volumes for the importance of study of third order structure functions for demys- 
tifying the two-dimensionalisation effect[8, 9] of the 3D turbulent fluid due to rapid rotation. 
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